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ABSTRACT  The  membrane-active  enzyme  phospholipase  D  (PLD)  catalyzes  the  hydrolysis  of  the  phosphodiester  bond  in 
phospholipids  and  plays  a  critical  role  in  cell  signaling.  This  catalytic  reaction  proceeds  on  lipid-water  interfaces  and  is  an 
example  of  heterogeneous  catalysis  in  biology.  Recently  we  showed  that  planar  lipid  bilayers,  a  previously  unexplored  model 
membrane  for  these  kinetic  studies,  can  be  used  for  monitoring  interfacial  catalytic  reactions  under  well-defined  experimental 
conditions  with  chemical  and  electrical  access  to  both  sides  of  the  lipid  membrane.  Employing  an  assay  that  relies  on  the 
conductance  of  the  pore-forming  peptide  gramicidin  A  to  monitor  PLD  activity,  the  work  presented  here  reveals  the  kinetics 
of  hydrolysis  of  long-chain  phosphatidylcholine  lipids  in  situ.  We  have  developed  an  extension  of  a  basic  kinetic  model  for  inter¬ 
facial  catalysis  that  includes  product  activation  and  substrate  depletion.  This  model  describes  the  kinetic  behavior  very  well  and 
reveals  two  kinetic  parameters,  the  specificity  constant  and  the  interfacial  quality  constant.  This  approach  results  in  a  simple  and 
general  model  to  account  for  product  accumulation  in  interfacial  enzyme  kinetics. 

INTRODUCTION 

Phospholipases  are  interfacial  enzymes  that  catalyze  the 
hydrolysis  of  ester  bonds  in  phospholipids.  These  enzymes 
play  an  important  role  in  lipid  metabolism,  cell  signaling, 
meiosis,  and  vesicle  trafficking  (1-5).  Due  to  the  amphi¬ 
philic  nature  of  their  phospholipid  substrates,  the  catalytic 
reaction  of  phospholipases  proceeds  on  membrane  inter¬ 
faces  and  depends  strongly  on  the  structure  and  properties 
of  these  interfaces  (1).  The  underlying  variables  that  govern 
catalysis  in  heterogeneous  environments,  such  as  the  phys¬ 
ical  structure  and  chemical  properties  of  lipid  interfaces, 
are  far  more  complex  than  those  encountered  by  soluble  en¬ 
zymes  in  homogeneous  solutions  (1-3, 6,7).  This,  in  turn, 
gives  rise  to  a  rich  kinetic  behavior  (8)  that  cannot  be 
described  with  simple  Michaelis-Menten  kinetic  analysis. 

Within  the  past  four  decades,  several  kinetic  models  have 
been  developed  to  describe  the  catalytic  reaction  of  interfa¬ 
cial  enzymes.  Verger  et  al.  pioneered  this  field  by  proposing 
the  first  and  simplest  kinetic  model  by  combining  the 
Michaelis-Menten  model  with  interfacial  activation  of 
enzymes  (9).  Since  then,  kinetic  models  have  been  proposed 
for  various  interfacial  structures  including  lipid  monolayers, 
liposomes,  and  micelles  (3,6,10,11).  Examples  include  the 
surface  dilution  kinetics  developed  by  Dennis  and  col¬ 
leagues,  which  described  catalysis  on  mixed  micelles 
(3,10,12),  and  the  scooting  and  hopping  modes  of  enzyme 
action  proposed  by  Berg  and  co-workers  to  describe  catal¬ 
ysis  on  liposomes  (6,11).  These  theoretical  models  have 
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been  frequently  employed  for  analyzing  the  kinetics  of 
phospholipases  on  model  membranes  (1 ).  Most  of  the  previ¬ 
ous  work  in  this  area  has,  however,  focused  on  kinetic 
analysis  of  phospholipase  A  (PLA)  (3,6,9,11-23).  More 
recently,  phospholipase  D  (PLD)  has  attracted  attention 
due  to  its  critical  role  in  cellular  processes  such  as  signaling, 
exocytosis,  and  migration  (24-32).  Only  a  limited  number 
of  studies  have  presented  any  analysis  of  the  interfacial 
kinetics  of  PLD  from  mammalian  (2,25,33)  and  plant 
(24,34,35)  cells,  and  most  of  these  studies  have  examined 
the  activity  of  PLD  on  mixed  micelles  and  employed  the 
surface  dilution  model  for  the  kinetic  analysis.  These 
models  typically  assume  that  the  products  of  the  enzyme  re¬ 
action  are  soluble  in  water  and  do  not  account  for  accumu¬ 
lation  of  phospholipid  products  with  long  acyl  chains. 

Here,  we  present,  to  our  knowledge,  the  first  quantitative 
kinetic  description  of  PLD  activity  on  planar  lipid  bilayers 
(36-40)  composed  of  long-chain  phospholipids  in  an 
attempt  to  mimic  the  physiological  conditions  of  long-chain 
phospholipid  substrates  and  products  that  are  associated 
with  cellular  membranes.  Starting  with  the  kinetic  model 
proposed  by  Verger  et  al.  (9)  for  short-chain  lipids,  we 
extend  this  model  to  account  for  the  interaction  between 
PLD  and  its  reaction  product,  phosphatidic  acid  (PA),  which 
is  a  long-chain  lipid  and  remains  in  the  membrane.  This 
analysis  also  demonstrates  that  a  recently  introduced  ion- 
channel-based  assay  (41)  that  reports  PLD-induced  changes 
of  the  ion  conductance  through  gramicidin  A  (gA)  pores 
(Pig.  1)  can  be  used  to  determine  the  kinetics  of  PLD-cata- 
lyzed  reactions.  It  should  be  noted  that  in  addition  to 
conductance,  PLD  activity  might  influence  other  aspects 
of  gA  channel  activity,  such  as  lifetime,  through  changes 
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FIGURE  1  Schematic  illustration  of  the  basic  concept  of  the  gA-based 
assay  employed  to  monitor  PLD  activity  (41).  Enzymatic  activity  of  PLD 
(green)  on  a  planar  lipid  bilayer  composed  of  electrically  neutral  PC  lipids 
(gray)  results  in  production  of  negatively  charged  PA  lipids  ( red );  the 
concomitant  accumulation  of  cations  (light  blue  spheres)  near  the  mem¬ 
brane  surface  leads  to  an  increase  in  single-channel  conductance,  y,  of  tran¬ 
sient  openings  of  gA  pores.  The  schematic  current- versus-time  traces  below 
the  cartoons  depict  the  PLD-induced  increase  in  the  ion  current  passing 
through  single  gA  channels. 


in  membrane  physical  properties  including  surface  charge, 
fluidity,  thickness,  and  curvature  (42-44).  Although  these 
aspects  of  gA  activity  might  provide  further  insight  into 
the  effect  of  PLD  binding  and  activity  on  membranes,  the 
assay  introduced  here  relies  solely  on  PLD-induced  changes 
in  membrane  surface  charge  that  can  be  monitored  through 
gA  conductance  and  does  not  consider  other  aspects  of  the 
gA  response  to  PLD  activity. 

MATERIALS  AND  METHODS 
Materials 

We  purchased  cesium  acetate,  cyclosporin  A,  and  1-octadecanethiol  from 
Sigma  Aldrich  (St.  Louis,  MO);  cesium  chloride  (CsCl)  from  International 
Biotechnologies  (New  Haven,  CT);  and  calcium  chloride  (CaCl2),  pentane, 
and  hexadecane  from  Fluka.  gA  was  purchased  as  gramicidin  D  from 
Sigma  Aldrich  and  purified  by  silica  chromatography  as  described 
previously  (45)  to  afford  a  final  purity  of  97%  of  gA.  We  purchased  the 
following  phospholipids  from  Avanti  Polar  Lipids  (Alabaster,  AL):  1,2-di- 
phytanoyl-5«-glycero-3-phosphocholine  (DiPhyPC),  1 ,2-diphytanoyl-sn- 
glycero- 3 -phosphate  (sodium  salt)  (DiPhyPA).  PLD  from  cabbage  (EC 
3. 1.4.4)  was  obtained  from  Sigma  Aldrich. 


Storage  and  final  concentration  of  PLD 

We  received  PLD  as  a  lyophilized  powder  and  dissolved  it  in  a  buffer  solu¬ 
tion  containing  10  mM  CsCl,  0.5  mM  CaCl2,  and  10  mM  cesium  acetate 
with  a  pH  of  5.5  (the  same  buffer  was  used  for  single-channel  recordings 
with  this  enzyme)  to  a  final  activity  of  2500  units  x  mL-1.  We  aliquoted 
and  stored  this  PLD  solution  at  — 80°C  until  usage.  According  to  Sigma 
Aldrich,  1  unit  of  PLD  liberates  1 .0  ^mol  of  choline  from  L-R-phosphati- 


dylcholine  (egg  yolk)  per  hour  at  pH  5.6  at  30° C.  The  specific  activity 
of  this  enzyme,  provided  by  Sigma  Aldrich,  was  >1670  units/mg  of 
protein.  Assuming  a  pure  enzyme  and  considering  a  molecular  weight  of 
~92,000  Da  (46),  a  concentration  of  1  unit  mL-1  corresponds  to  a  PLD 
concentration  of  ~6.5  nM. 


Formation  of  planar  lipid  bilayers 

We  formed  planar  lipid  bilayers  with  the  folding  technique  (38),  as 
described  previously  (41).  Briefly,  two  compartments  (with  3-  or  4-mL 
capacity)  of  a  custom-made  Teflon  chamber  were  separated  by  a  thin  Teflon 
film  containing  a  single  pore  with  a  diameter  of  -100  pm.  To  facilitate 
bilayer  formation,  we  pretreated  the  area  around  the  pore  with  2  pL  of 
5%  (v/v)  hexadecane  in  pentane.  Upon  addition  of  a  fraction  of  electrolyte 
solution  (~  1  mL)  to  each  compartment,  we  spread  3-5  pL  of  a  solution  of 
25  mg  mL- 1  DiPhyPC  in  pentane  at  the  air- water  interface.  Then  we  raised 
the  liquid  level  by  adding  the  rest  of  the  electrolyte  solution  (2  or  3  mL)  to 
each  compartment,  which  resulted  in  bilayer  formation  as  described  origi¬ 
nally  by  Montal  and  Mueller  (47).  We  repeated  the  cycle  of  raising  and 
lowering  the  liquid  levels  until  we  obtained  a  bilayer  with  a  minimum 
capacitance  of  70  pF. 


Single-channel  recordings 

After  formation  of  a  stable  lipid  bilayer,  we  sequentially  added  small 
volumes  (0.2  pL)  of  a  solution  of  10  ng  mL-1  gA  in  isopropanol  (Acros 
Organics,  Geel,  Belgium)  to  both  compartments  of  the  bilayer  setup  until 
one  to  six  active  gA  channels  could  be  observed  in  the  bilayer  at  the 
same  time.  After  each  addition  of  gA,  we  mixed  the  bilayer  chambers  by 
stirring  the  solutions  in  both  compartments  for  at  least  3  min  with  a  stir 
bar  (Sun  Stir,  Warner  Instruments,  Hamden,  CT).  In  these  experiments, 
the  final  concentration  of  gA  in  the  bilayer  chamber  was  in  the  range 
0. 1-2.0  pM.  This  range  of  gA  concentration  corresponds  to  a  total  amount 
of  0.3-8  fmol  of  gA  in  the  chamber  with  a  volume  of  3^1  mL;  the  total 
amount  of  lipid  in  each  compartment  was  -150  nmol,  leading  to  a  lipid/ 
gA  ratio  of  -108.  To  measure  the  single-channel  conductance  of  gA  pores 
(48,49),  we  recorded  current- versus-time  traces  while  applying  different 
voltages  in  the  range  of  ±100  mV  across  the  planar  lipid  bilayers.  We  per¬ 
formed  these  single-channel  recordings  in  voltage-clamp  mode  using  Ag / 
AgCl  pellet  electrodes  (Warner  Instruments)  in  both  compartments  of  the 
bilayer  setup.  Data  acquisition  and  storage  were  carried  out  using  custom 
software  written  in  Lab  view  by  Dr.  J.  D.  Uram  and  Dr.  D.  J.  Estes  in  com¬ 
bination  with  a  Geneclamp  500  amplifier  from  Axon  Instruments  (Union 
City,  CA;  set  to  a  gain  of  100  mV  pA_  and  a  filter  cutoff  frequency  of 
2  kHz).  The  data  acquisition  board  (National  Instruments,  Austin,  TX) 
that  was  connected  to  the  amplifier  was  set  to  a  sampling  frequency  of 
15  kHz.  All  recordings  were  carried  out  at  a  temperature  of  ~22°C.  Record¬ 
ings  could  be  routinely  performed  for  15-20  min;  due  to  limited  stability  of 
the  bilayer,  however,  recording  for  longer  times  became  increasingly  more 
challenging  as  the  experiment  progressed. 

We  analyzed  single-channel  current  traces  by  computing  histograms  of 
the  currents  from  the  original  current-versus-time  traces  with  Clampfit 
9.2  software  (Axon  Instruments)  (41).  From  these  histograms,  we  extracted 
the  mean  current  amplitude  of  gA  channel  opening  and  closing  events.  All 
conductance  values  were  obtained  from  the  slopes  of  current-amplitude- 
versus-voltage  (I-V)  curves.  Fig.  1  demonstrates  the  experimental  setup 
used  for  this  study. 


Preparation  of  small  liposomes 

We  prepared  small  unilamellar  liposomes  of  DiPhyPC  lipids  by  tip  sonicat- 
ion  as  described  previously  (50,51).  Briefly,  we  deposited  a  small  droplet 
of  the  lipid  stock  solution  (25  mg  mL-  DiPhyPC  in  pentane)  into  a 
clean  round-bottom  flask  with  a  volume  of  5  mL  and  employed  a  rotatory 
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evaporator  to  form  a  thin  lipid  film  on  the  wall  of  the  flask.  To  remove  re¬ 
sidual  traces  of  pentane,  we  desiccated  the  lipid  film  under  vacuum  (—740 
torr)  for  at  least  1  h.  We  hydrated  the  lipid  film  by  adding  an  aqueous  solu¬ 
tion  of  150  mM  KC1  to  reach  a  final  lipid  concentration  of  6  mM.  We 
formed  small  liposomes  by  tip-sonicating  the  resulting  solution  using  a 
Branson  Sonifier  150  (Branson  Ultrasonics,  Danbury,  CT)  for  10-12  min 
(with  ~5  W  output  energy).  During  tip  sonication,  the  flask  was  immersed 
in  an  ice  bath  to  prevent  excessive  heating  of  the  solution. 


Surface  plasmon  resonance  measurements 

We  prepared  glass  substrates  to  use  for  surface  plasmon  resonance  (SPR) 
measurements  from  SF10  glass  slides  (n  —  1.7)  (Accurion,  Goettingen, 
Germany).  We  coated  one  side  of  these  slides  with  a  2-nm-thick  layer  of 
titanium  followed  by  a  50-nm-thick  layer  of  gold.  Immediately  after  coating 
the  slides,  we  immersed  them  overnight  in  a  solution  containing  5  mM  oc- 
tadecanethiol  to  generate  a  self-assembled  monolayer  (SAM)  of  octadecane 
on  the  gold  surface  of  the  slide.  We  used  a  Nanofilm  EP3-SPR  instrument 
(Accurion)  with  an  SPR/total-intemal-reflection  fluid-cell  accessory  for 
the  SPR  measurements.  We  mounted  the  SF10  glass  slide  in  the  fluid  cell 
with  a  60°  SF10  prism  and  index  matching  fluid  (n  —  1 .7)  between  the  slide 
and  the  prism.  We  set  the  instrument  to  direct  polarized  light  (A  =  532  nm) 
through  the  prism  onto  the  gold-coated  chip  at  an  incident  angle  of  60°. 
Before  each  experiment,  we  degassed  all  buffer  solutions  to  prevent  air  bub¬ 
bles  from  forming  in  the  fluid  cell.  We  flowed  solutions  through  the  fluid 
cell  and  over  the  octadecane  monolayer,  proceeding  according  to  the 
following  order:  1),  Filling  the  cell  with  reaction  buffer  by  flowing  a  solu¬ 
tion  of  PLD  buffer  (10  mM  CsCl,  0.5  mM  CaCl2,  and  10  mM  cesium 
acetate,  pH  5.5)  at  a  flow  rate  of  100  fih  min-1  for  5  min;  2)  generating 
a  self-assembled  lipid  monolayer  on  top  of  the  octadecane  SAM  by  intro¬ 
ducing  a  solution  of  liposomes  containing  2  mM  DiPhyPC  and  150  mM 
KC1  at  200  fiL  min-1  for  15  min;  3)  removal  of  excess  liposomes  by  flowing 
a  solution  of  PLD  buffer  at  100  /iL  min-1  for  20  min;  4),  measuring  the  on 
rate  of  PLD  binding  to  the  lipid  layer  by  introducing  PLD  buffer  that  con¬ 
tained  the  PLD  (at  a  concentration  of  20,  41,  or  1 10  nM)  at  100  /xL  min-1 
for  20-30  min;  and  5),  measuring  the  off  rate  of  PLD  from  the  lipid  layer  by 
introducing  PLD  buffer  (without  the  enzyme)  at  100  /iL  min-1  for  30  min. 
Fig.  SI  in  the  Supporting  Material  illustrates  a  representative  SPR  signal 
during  an  experiment.  During  each  experiment,  we  monitored  the  change 
in  the  Delta  signal,  A D,  and  the  Psi  signal,  Ai|/,  from  the  instrument  using 
the  EP3  View  Software  (Accurion).  From  these  signals,  we  calculated 
the  exponential  rate  constant,  t,  for  each  concentration  of  PLD  from: 
A D(t)  =  A Dmax  x  (1  -  exp(r/r))  and  A \f/  ( t )  =  A \J/max  x  (1  -  exp(f/r)) 
(see  Fig.  S2  A).  The  ratio  1/r  is  linearly  dependent  on  the  bulk  concentration 
of  PLD  as  described  by  the  equation  1/r  =  kon  [PLD]  +  kof[  (52).  Therefore, 
we  plotted  (1/r)  as  a  function  of  [PLD]  and  obtained  the  on  rate  ( kon )  and  off 
rate  (£0ff)  from  the  best  linear  fit  to  the  data  in  this  plot  (see  Fig.  S2  B ).  The 
slope  of  the  linear  fit  to  the  data  corresponded  to  kon  and  the  intercept  of  the 
linear  fit  with  the  y  axis  corresponded  to  kQff.  We  used  these  values  to  calcu¬ 
late  Kd  —  koff/kon.  Note,  kon  and  koff  in  these  SPR  experiments  correspond  to 
kp  and  kd,  respectively,  in  the  kinetic  analysis. 


Phosphatidylcholine  Choline  PhosphatidicAcid 

(PC)  (PA) 


FIGURE  2  (A)  Schematic  of  the  enzymatic  reaction  catalyzed  by  PLD. 
( B )  Illustration  of  an  extended  version  of  Verger’s  Michaelis-Menten  model 
(9)  to  describe  the  interfacial  hydrolysis  of  lipids.  Illustrated  steps  include 
reversible  adsorption  of  the  soluble  enzyme  (E)  to  an  interfacial  form  (E*) 
followed  by  a  2D  Michaelis-Menten  reaction  where  the  enzyme  binds  the 
interfacial  substrate  (S*)  to  form  the  Michaelis  complex  (E*S*)  and  then 
releases  the  product(s),  which  may  remain  in  the  membrane  (Pi*)  and/or 
be  released  back  into  solution  (P2).  Although  the  reaction  products  of 
membrane-active  enzymes  are  often  long-chain  lipids  and  therefore  not 
water-soluble.  Verger’s  model  is  best  suited  for  short-chain  lipids  whose 
hydrolysis  products  can  be  treated  as  soluble  species  that  diffuse  away 
from  the  lipid  membrane  into  the  aqueous  bulk  solution  (1).  The  cartoon 
shown  here  is  an  extended  version  of  Verger’s  original  cartoon  and  depicts 
the  possible  mechanisms  {gray  arrows )  by  which  product  accumulation  at 
the  interface  might  modify  the  kinetics  of  individual  steps  along  the  inter¬ 
facial  catalytic  reaction. 


leases  the  product(s),  which  may  remain  in  the  membrane 
(P*)  and/or  be  released  back  into  solution  (P).  By  invoking 
the  usual  steady-state  assumptions,  Verger  et  al  arrived  at 
the  following  equation  for  the  velocity  (9) 


d  rv 

dt 


*2  rs* 

- ACl 

(^M  +  rS*)  y 


0) 


RESULTS  AND  DISCUSSION 
Kinetic  model 

Verger  et  al.  proposed  the  first  adaptation  of  the  Michaelis- 
Menten  model  to  describe  the  interfacial  kinetics  of  PLA 
(9).  Fig.  2  illustrates  the  general  scheme:  the  basic  steps 
involve  reversible  adsorption  of  the  soluble  enzyme  (E)  to 
an  interfacial  form  (E*  )  followed  by  a  2D  Michaelis-Menten 
reaction  where  the  enzyme  binds  the  interfacial  substrate 
(S*)  to  form  the  Michaelis  complex  (E*S*)  and  then  re- 


where  KM*  =  (k2  +  k  _ ,  )/k ,  is  the  interfacial  Michaelis- 
Menten  constant,  rP*  and  Ts*  (mol  m  2)  are  the  surface 
concentrations  of  the  product  and  the  substrate,  respectively, 
CEo  (mol  in  3)  is  the  total  enzyme  concentration,  A  (m2)  is 
the  total  water/lipid  interfacial  area  (in  our  setup  -2.25  x 
I  IF  1  m2),  and  V  (m3)  is  the  volume  of  the  electrolyte  solu¬ 
tion  (in  our  setup  3  or  4  x  10-6  m3). 

This  model  assumes  that  the  product  P*  diffuses  quickly 
and  does  not  affect  subsequent  reactions.  This  assumption 
may  be  valid  for  water-soluble  short-chain  phospholipids 
but  is  likely  not  true  for  natural  long-chain  phospholipids. 
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which  remain  in  the  bilayer  (13).  In  fact,  accumulation  of  re¬ 
action  products  at  the  interface  can  dilute  the  substrate  at  the 
surface,  inhibit  the  enzymatic  reaction,  or  even  alter  the 
binding  of  enzyme  to  the  interface  and  thus  affect  the  hydro¬ 
lysis  reaction  (1,9).  The  gray  arrows  in  Fig.  2  illustrate  how 
these  effects  could  modify  the  kinetics  of  various  steps  in 
the  enzymatic  reaction. 

In  the  case  of  PLD,  we  find  exactly  these  effects.  First,  as 
the  product  PA  accumulates  in  the  membrane,  we  observe 
an  increase  in  PLD  activity,  suggesting  that  the  product  is 
affecting  the  equilibrium  between  E  and  E*,  represented 
by  kjk p  in  Eq.  1.  A  cation-mediated  interaction  between 
PA  and  PLD  has  been  documented  (53),  and  we  account 
for  this  effect  by  modeling  the  adsorption-desorption  of 
the  enzyme  as  a  bimolecular  association  process, 

kp 

P*  +E  ^  E*, 

ki 

where  E*  represents  the  membrane-associated  enzyme  bond 
to  PA  lipids.  As  a  result,  the  equilibrium  between  E  and  E* 
in  the  denominator  of  Eq.  1  changes  from  kd/kp  to  kd/(kprP*) 
(see  Appendix  for  a  derivation  of  a  more  complete  model). 
The  second  contribution  that  we  need  to  account  for  is  sub¬ 
strate  depletion  and  dilution  due  to  product  accumulation. 
The  lipid  density  of  the  bilayer  is  fixed  (ro)  and  is  simply 
the  sum  of  the  substrate  and  product  densities,  assuming 
that  the  area/lipid  is  approximately  the  same  for  substrate 
and  product  (54) 

T0  =  IV  +  rP*. 


One  further  step  can  be  taken  to  simplify  this  model,  since 
for  the  bilayer  system  used  in  this  work  the  A/V  ratio  is  small 
and  the  second  term  in  the  denominator  of  Eq.  1  can  be 
ignored.  Taken  together,  we  arrive  at  an  extended  descrip¬ 
tion  for  the  velocity: 


d  TP. 
dt 


— (r'o  -  rP.)rP*cEo  —  Gm(r0  —  rP*)rP»cE(), 

r*  kd 
Kuk 

Kp 

(2) 


where  we  have  introduced  the  interfacial  quality  term  Qm 
(9).  Equation  2  gives  the  rate  of  product  formation,  and  inte¬ 
grating  this  expression  gives  the  product  concentration  as  a 
function  of  time, 


rP.  = 


r 


o 


1  +  a  exp(  -  0mroCEo?)’ 


(3) 


or,  written  as  a  dimensionless  mole  fraction. 


-XpA  — 


1 

1  +  a  exp(  -  GmroCEof)’ 


(4) 


where  a  is  an  integration  constant. 


Quantifying  PLD  activity 

We  monitored  the  enzymatic  activity  of  PLD  using  the  pre¬ 
viously  described  gA-based  assay  (41).  Briefly,  we  formed 
planar  lipid  bilayers  composed  of  long-chain  PC  lipids 
and  added  the  appropriate  concentration  of  gA  (0.1— 2.0 
pM)  to  the  bulk  electrolyte  solution  to  generate  only  a 
few,  transiently  open  gA  pores  at  any  given  time  in  these 
bilayers.  We  added  PLD  to  the  electrolyte  solutions  on 
both  sides  of  the  lipid  bilayer  and  monitored  the  changes 
in  the  single-channel  conductance  of  gA  pores,  y,  over 
time.  As  expected,  the  enzymatic  activity  of  PLD  on  the 
planar  lipid  bilayer  led  to  a  time-dependent  increase  in 
the  single-channel  conductance  of  gA.  This  increase  is  the 
result  of  PLD-mediated  formation  of  negatively  charged 
PA  lipids  in  the  membrane,  which  leads  to  accumulation 
of  cations  close  to  the  entrance  of  gramicidin  pores  in  the 
membrane.  Since  monovalent  cations  are  the  charge  carriers 
that  translocate  through  gA  pores,  this  local  accumulation  of 
ions  leads  to  a  measurable  increase  in  y.  Fig.  3  demonstrates 
this  increase  in  y  as  a  function  of  time  after  the  addition  of 
PLD  for  three  different  enzyme  concentrations. 

To  quantify  the  enzymatic  activity  of  PLD,  we  prepared  a 
calibration  curve  of  y  as  a  function  of  the  mole  fraction  of 
PA  lipids  as  previously  described  (41).  From  this  calibration 
curve,  which  was  done  in  the  same  buffer  as  the  experiments 


A 


C 


FIGURE  3  ( A  and  B )  Changes  in  the  single¬ 

channel  conductance  of  gA,  y,  after  the  addition 
of  PLD.  (A  and  B )  Current  versus  time  recordings 
in  the  bilayer  containing  gA  before  addition  of 
PLD  (A)  and  13  min  after  addition  of  20.3  nM 
PLD  ( B ).  (C)  Time-dependent  change  in  gA 
conductance  after  the  addition  of  15.2  ( black 
ovals),  20.3  ( red  circles),  and  33.3  nM  (, green 
diamonds)  concentrations  of  PLD  to  both  sides  of 
the  bilayer.  Error  bars  represent  the  mean  ±  SE 
(N  >  3). 
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with  PLD  (see  Supporting  Material),  we  obtained  Eq.  5, 
which  describes  the  relationship  between  y  (in  units  of 
pS  =  10-12  fi-1)  and  the  mole  fraction  of  PA  in  the 
bilayer  (XPA) 


y  =  3.50  pS  + 


12.83  pS  x  Apa 
0.23  +  XpA 


(5) 


Using  Eq.  5,  we  converted  the  time-dependent  changes 
in  y  upon  addition  of  PLD  to  time-dependent  changes  of 
the  mole  fraction  of  PA  lipids,  XPA,  in  the  membrane. 
Fig.  4  A  shows  plots  of  XPA  as  a  function  of  time  for  PLD 
concentrations  of  Ce0  =  15.2,  20.3,  27.1,  and  33.3  nM.  In 
this  figure,  data  are  fit  using  the  novel  kinetic  model  in 
Eq.  4.  Fig.  4  B  shows  the  same  data  after  rescaling  the 
time  axis  to  (Ce0  x  time).  According  to  Eq.  4,  the  kinetic 
model  developed  here  predicts  that  this  rescaling  should 
result  in  all  of  the  data  collapsing  to  a  single  curve,  and 
indeed  we  find  no  statistical  difference  between  the  data 
and  this  theoretical  prediction  (Kolmogorov-Smirnov  test 
D  =  0.1071,  p  =  0.9971,  N  =  28). 

Two  kinetic  parameters,  the  specificity  constant,  k2/KM*, 
and  the  interfacial  quality  constant,  Qm  =  k2l  (Km*  x  (kj 
A'p)),  are  commonly  used  to  describe  enzymatic  reaction 
rates  of  phospholipases  quantitatively.  Qm  has  slightly 
different  units  in  our  model  since  PLD  association  with 
the  membrane  is  now  lipid  dependent,  and  from  curve  fitting 
in  Fig.  4,  we  obtained  a  value  of  14.5  ^M-1  s-1  for  Qm  ro  of 
the  examined  PLD  activity  on  long-chain  PC  lipids.  To 
estimate  the  specificity  constant,  we  first  employed  SPR  to 
determine  the  binding  affinity  of  PLD  to  PC  bilayers  and 
obtained  the  dissociation  constant  for  this  binding  event  in 
units  of  bulk  concentrations,  Kd  =  150  nM.  This  number, 
based  on  the  definition  of  Kd,  represents  the  enzyme  concen¬ 
tration  (Ceo)  at  which  half  of  the  available  binding  sites  for 
enzyme  are  occupied  by  bound  enzymes.  We  thus  assume 
that  at  the  enzyme  concentration  of  CEo  =150  nM,  half 
of  the  maximally  available  surface  area  of  the  lipid  bilayer 
is  occupied  by  bound  PLD.  To  determine  the  surface  con¬ 
centration  of  the  enzyme,  rE*,  when  Ce0  =  Kd,  we  approx¬ 
imated  the  size  of  individual  PLD  molecules  as  described 


below,  and  thus,  we  could  determine  the  equilibrium  disso¬ 
ciation  constant,  kd/(kp  ro)  =  CEo/rE*. 

We  approximated  the  size  of  individual  PLD  proteins 
employed  in  the  assay  presented  here  based  on  the  molecular 
weight  of  -92  kDa  (55)  and  a  comparison  with  another  pro¬ 
tein  with  the  same  molecular  weight  (a  92-kDa  fragment  of 
yeast  DNA  topoisomerase  II)  (56),  whose  molecular 
volume  of  720  nm3  is  known.  Assuming  a  spherical  shape, 
each  PLD  molecule  would  cover  an  area  of  -97  nm2  on  the 
bilayer.  If  we  consider  a  square  unit  area  of  1  pm2  and  a  close 
hexagonal  packing  (corresponding  to  a  coverage  of  90.7%) 
(57)  of  the  PLD  enzyme  at  the  interface,  the  maximum  sur¬ 
face  coverage  would  be  -9.4  x  103  enzymes  pm~  2  at  the 
interface.  This  coverage  corresponds  to  the  densest  possible 
packing  of  PLD  in  the  absence  of  any  steric  hindrance  or  ther¬ 
mal  motion.  If  we  consider  less  ideal  packing,  where  each 
PLD  molecule  would  occupy  twice  its  actual  projected  area 
on  the  membrane,  this  would  reduce  the  maximal  possible 
coverage  to  -45%  or  -4.7  x  103  enzymes  /inT  2  at  the  mem¬ 
brane  surface.  Thus,  we  assume  that  maximal  PLD  coverage 
on  the  membrane  would  be  within  the  range  45-90.7%.  At 
half-maximal  binding,  therefore, -2.3  x  103-4.7  x  103  en¬ 
zymes  pmT2  could  be  bound  at  the  interface,  corresponding 
to  an  approximate  surface  concentration  of  TE*  =  3.8  - 
7.8  nmol  enzyme  m-2.  Consequently,  the  constant  CEo/rE* 
and,  hence,  kd/(k„  ro)  have  an  estimated  value  of  19.0  x 
103  -  39.5  x  103  m  1 .  Based  on  this  range  of  values  for 
kd/(kp  ro),  we  calculated  a  range  of  values  from  our  experi¬ 
mentally  determined  value  of  Qm  ro  for  the  specificity 
constant,  k2/KM*,  of  0.77  -  1.6  m2  nmol-1  s-1  for  the 
PLD-catalyzed  hydrolysis  of  long-chain  PC  lipids. 

PLD-catalyzed  hydrolysis  of  long-chain  PC  lipids  results 
in  production  of  long-chain  PA  lipids,  which  remain  in  the 
membrane.  Fig.  4  clearly  illustrates  that  PLD  activity  on 
the  PC  bilayer  starts  with  a  relatively  slow  phase,  and  as 
PA  accumulates  in  the  membrane,  PLD  catalytic  activity 
accelerates.  Accumulation  of  PA  in  the  presence  of  Ca2+, 
which  is  required  for  maximal  activity  of  cabbage  PLD, 
leads  to  the  formation  of  Ca2+/PA  domains  in  the  membrane 
(41,53).  Kuppe  et  al.  have  previously  shown  that  these 


Time  (min) 


Scaled  Time  (time  x  [PLD]) 


FIGURE  4  Time  dependence  of  PA  production 
for  different  concentrations  of  PLD.  (A)  The 
mole  fraction  of  PA  (XpA)  for  PLD  concentrations 
of  15.2,  20.3,  27.1,  and  33.3  nM.  ( B )  The  same  data 
as  in  A  after  rescaling  the  time  by  the  PLD  concen¬ 
tration.  As  predicted  by  Eq.  4,  the  data  points 
collapse  to  a  single  curve.  The  magenta  curve  is 
the  full  model  from  Eq.  8  and  is  indistinguishable 
from  the  simplified  model  given  by  Eq.  4  ( black 
curve).  The  error  bars  in  both  plots  represent  the 
mean  ±  SD  (N  >  3)  and  the  curves  are  best  fits 
to  the  data  using  Eq.  4. 
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domains  facilitate  the  binding  of  PLD  to  the  membrane  and 
hence  lead  to  an  increase  in  the  hydrolysis  rate  of  this 
enzyme  (53).  Ahyayauch  et  al.  suggested  that  Ca2+/PA 
complexes  bind  to  an  allosteric  site  of  PLD  and  facilitate 
its  association  with  the  membrane.  Moreover,  the  activity 
of  phospholipases  is  strongly  influenced  by  physical  proper¬ 
ties  of  the  interface  such  as  membrane  curvature  (58).  PA  is 
a  cone-shaped  lipid  and  is  known  to  induce  negative  curva¬ 
tures  in  bilayers  that  could  contribute  to  modulation  of  PLD 
binding  and  activity. 

A  similar  initial  slow  phase,  typically  referred  to  as  a 
lag  phase,  has  previously  been  reported  for  other  phospholi¬ 
pase  enzymes  (59).  As  PLD  activity  further  proceeds  on 
the  membrane,  however,  accumulation  of  product  leads  to 
depletion  of  substrate  lipids  in  the  membrane  and  thus  slows 
the  PLD  catalytic  rate.  The  kinetic  model  presented  here 
considers  these  two  effects  of  PA  accumulation  and,  as 
demonstrated  in  Fig.  4,  describes  the  behavior  of  this 
enzyme  extremely  well. 

The  specificity  constant  and  interfacial  quality  constant, 
the  two  parameters  that  we  determined  in  this  work,  are 
commonly  employed  to  describe  the  kinetics  of  interfacial 
enzymatic  reactions  (1,9).  The  specificity  constant  is  a 
measure  of  the  catalytic  efficiency  of  the  enzyme,  and  the 
interfacial  quality  constant  accounts  for  the  effect  of 
physiochemical  parameters  of  the  interface,  such  as  the 
orientation  and  conformation  of  the  molecules,  the  molecu¬ 
lar  and  charge  densities,  the  water  structure,  and  the  vis¬ 
cosity,  on  the  activity  of  the  enzyme  (1).  The  value  of 
14.5  /iM_ 1  s  1  for  the  interfacial  quality  constant  obtained 
for  the  enzymatic  reaction  of  PLD  on  the  long-chain  PC 
lipid  examined  here  appears  reasonable  for  this  enzyme 
based  on  a  previously  reported  value  of  3  /iM"  1  s-1  for 
the  hydrolysis  of  a  medium-chain  length  lipid  by  a  different 
lipase  enzyme,  PLA,  on  a  lipid  monolayer  (60). 


chemical  access  to  both  sides  of  the  lipid  bilayer.  In  addi¬ 
tion,  planar  lipid  bilayers  can  be  prepared  from  a  variety 
of  lipids  (62-64).  Likewise,  the  theoretical  description  of 
PLD  activity  introduced  here  will  contribute  to  the  under¬ 
standing  of  the  collective  activity  of  systems  of  lipolytic 
enzymes  and  their  function  in  cell  signaling. 


APPENDIX 

The  kinetic  model  we  present  that  results  in  Eq.  4  is  based  on  the  assump¬ 
tion  that  PA  in  the  membrane  attracts  more  PLD  to  the  surface.  We  account 
for  this  effect  by  including  a  bimolecular  association  reaction  where  PLD 
(E)  binds  to  PA  (P*)  to  produce  a  membrane-associated  form  of  PLD 
(E*).  This  is  a  simple  model,  but  technically  (and  mathematically)  it  has 
problems  if  we  start  with  a  pure  PC  membrane,  since  without  any  PA  no 
enzyme  will  ever  associate  with  the  membrane.  A  more  correct  model 
would  adopt  a  form  for  the  penetration  rate  constant  like 

*„  =  *p(i  +  /nv),  (6) 


where  we  now  have  some  basal  rate  of  PLD  association  to  the  membrane 
(kp)  that  naturally  increases  as  more  PA  accumulates  (kpfiYp*).  When  we 
adopt  this  form,  we  arrive  at  a  more  complex  form  for  the  velocity. 


dTp* 

~dT 


Qm(T0  -  I»(l  +  prP>)CEo. 


(7) 


If  we  integrate  Eq.  7  and  impose  the  condition  that  we  initially  have  no 
PA  in  the  membrane  (at  t  —  0),  we  arrive  at  a  slightly  more  complex  formula 
for  the  product  concentration  as  a  function  of  time, 

P  =  r0(exP[(i  +  ffr0)emcy]  -  l) 

P*  exp[(l  +  pr0)QmcEj]  +  jSro  ' 

We  fit  our  scaled  data  using  this  more  complete  model  and  get  the  plot 
shown  in  Fig.  4  B,  where  the  black  line  is  the  original  fit  from  Eq.  4  and 
the  magenta  curve  is  the  complete  model  based  on  Eq.  8.  The  difference 
between  the  simple  and  full  model  is  obviously  very  minimal  and  only 
occurs  for  small  values  of  Xpa-  Indeed,  the  best  fit  to  the  data  gives  a  value 
of  jd  =  178.7,  and  in  the  limit  of  (3  »  1,  we  can  ignore  the  1  in  the  numer¬ 
ator  and  Eq.  8  reduces  to  the  same  form  as  Eq.  3. 


CONCLUSIONS 

This  study  demonstrates  that  the  gA-based  assay  employed 
here  can  be  used  for  monitoring  the  kinetics  of  catalytic 
reactions  of  membrane-active  enzymes  on  natural  lipid 
substrates  whose  net  electrical  charge  changes  as  a  result 
of  enzyme  activity.  Further,  the  high-quality  data  arising 
from  this  assay  allowed  us  to  develop  a  novel  and  quantita¬ 
tive  model  that  describes  the  activity  of  interfacial  enzymes 
on  physiologically  relevant  long-chain  phospholipids  and 
accounts  for  the  accumulation  of  product  lipids  and  dilution 
of  substrate  lipids  in  the  membrane.  Based  on  this  analysis, 
we  were  able  to  determine  the  specificity  constant  and  inter¬ 
facial  quality  constant  for  the  examined  reaction.  Planar 
lipid  bilayers  provide  an  attractive  platform  for  such  studies, 
since  they  share  many  of  the  advantages  of  lipid  monolayer 
systems,  as  compared  to  the  bulk  systems  such  as  micelles 
and  liposomes  (61),  and  mimic  the  lipid  bilayer  structure 
of  natural  cell  membranes  while  providing  electrical  and 
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